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We investigate color superconductivity and chiral symmetry restoration at finite temper¬ 
ature and baryon density in the frame of standard two flavor Nambu-Jona-Lasinio model. 
We derive the diquark mass in RPA, discuss its constraint on the coupling constant in 
the diquark channel, and And a strong competition between the two phase transitions 
when the coupling constant is large enough. 
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1. Introduction 


In the ideal case at asymptotically high baryon density, the color superconductiv¬ 
ity with two massless flavors and the color-flavor-locking (CFL) phase with three 
degenerated massless quarks have been widely discussed from first principle QCD 
calculation^. For physical applications we are more interested in the moderate 
baryon density region which may be related to the neutron stars and, in very opti¬ 
mistic cases, even to heavy-ion collisions. Usually, four-fermion interaction models 
such as the instanton, as well as the Nambu-Jona-Lasinio (NJL) model, are used in 
the study of color superconductivity phase transition. It was found that the color 
superconductivity gap can be of the order of 100 Me^f^, which is two orders larger 
than early pertu rbative QCD estimatiorPl. The electric and color charge neutral¬ 
ity conditioJ^I^ leads to a new phase of color smerconductivity, the gapless color 
superconductivitjElor the breached pairing phasJH The most probable temperature 
for this new phase is finite but not zer<M 

It is generally accepted that the NJL model^ applied to quarkJl^ offers a simple 
but effective scheme to study spontaneous chiral symmetry breaking in the vacuum, 
chiral restoration at finite te mperature and density, and spontaneous color sym¬ 
metry breaking at high density^l^lsLyj From the study of chira l ph ase transition 
without considering color superconductivity, it is well-known thalM^ the mean field 
approximation to quarks and random phase approximation (RPA) to mesons can 
describe well the thermodynamics of the system, especially the massless Goldstone 
mode in the spontaneous symmetry breaking phase. In this letter I will determine 
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the diquark mass in RPA, discuss its constraint on the coupling constant in the 
diquark channel, and investigate the phase structure of color superconductivity and 
chiral restoration and its dependence on the coupling constant. 


2. Quarks in Mean Field Approximation 


The NJL Lagrange density is defined as 


£ =-m + ^7o)V’ (1) 

where Gs and Gd are coupling constants in color signet channel and color anti¬ 
triplet channel, respectively, and m is the current quark mass. 

The quark-antiquark and diquark condensates which are order parameters of 
chiral and color superconductivity phase transitions are defined as 

a = , A = A3 = ^ ^ ( 2 ) 


where it has been regarded that only the first two colors participate in the con¬ 
densate, while the third one does not. The condensates are assumed real but are 
otherwise as yet unspecified. 

A usual way to treat systematically the spontaneous color and chiral symme¬ 
try brea king at finite temperature and density is to introduce the Nambu-Gorkov 
spac^i^ with color and flavor degrees of freedom. After performing the Fourier 
transformation, the quark propagator in momentum space at mean field level is 
diagonal in the 12 dimensional Nambu-Gorkov space, 


S{k) = 


(SAik) \ 

SB{k) 

Scik) 

Snik) 

SE{k) 

V SF{k)J 


( 3 ) 


where S[{k) with I = A, B, C, D, E, F are 2x2 matrices. 


Si = 



( 4 ) 


In general case, the quark chemical potential fXia depends on color and flavor degrees 
of freedom, the elements of the 6 matrices Sj are totally different. However, if we 
consider only the baryon chemical potential one has ^ = ^_b/ 3, and some of 
the elements are identical , 


G + _ /^+ _ /^+ _ /^+ 

A ~ ~ ~ ’ 

_ _^+ _ _^+ 
~ 


G + _ /^+ 

E — ^ F 


G-E = G-p 


G-^ = G-^ = Gc = G-^, 


, — ^23 , 

“E ~ ~ “E = “E ~ ^ ’ 


I 
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( 5 ) 
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G+{k) 

G^{k) 

si(fc) 

^Aik) 

GUk) 

G^ik) 


/Cq 


^0 — -r 
{ko - Ef. ){ko + Ef. ) 


(fco ■ 
(fco ■ 


2iGD^ ^ 2iGr)^ ^ 

■E^Kko + Ej:)^ + (ko-El)iko + E+)^ 

2iGD^ 2iGDA 

- E^){ko + E^) + {ko - E+){ko - 

1 1 




1a 1a 

- -——A+7o + —A_7o 

ko ~ E}. + /r fco + Afc + /X 


1 ^ ^ 1 
fco - Efe - /X fco + Efc - /X 


■A-7o , 


with the quark energies E^ = {Ek ± /x)^ + (2 G'dA)^ , Ek = y^lkp + Mg, 
effective quark mass Mq = m — 2Gs<J , and the energy projectors A± = 

1 A ± 7o(7-k+MQ) \ 

2 y Ek ) ' 

While the first two colors are involved in the quark-antiquark and diquark con- 
densates, and their dispersion relation are related to a and A, the third color 
participates in the chiral condensate only, the non-diagonal elements of the last two 
matrices Se and Sp are zero, and its energy Ek is associated with cr only. 

From the definition of the two condensates ©, it is easy to obtain their relation 
with the elements of the quark propagator 


" = "«/tS; w (2/<En -1) 


-b 


(27r)3F;fcL E, 

{2fiE+) - 1) + f{Ek - m) + f{Ek + m) - 1 


( 6 ) 


A 


Ek + /i 

Et 

' + / (Sf {-k WEt) - 1) 


with the Fermi-Dirac distribution function f{x) = + 1). 


= 0 , 


3. Coupling Constant Gd in Diquark Channel 

The diquark polarization function can be represented in terms of the matrix ele¬ 
ments of the quark propagator, 

- illDik) = -^y [h5iG\{p - k)ij5iG^ip)] ■ (7) 

In the mean field approximation to quarks and RPA to mesons and diquarks, 
the temperature and baryon chemical potential dependence of the diquark mass 
Md{T,p) is controlled by the pole equation 

l-2GDnz7(fc2 = M|„O) = 0 . 


( 8 ) 
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We consider now the physical constraints in the vacuum on the coupling constant 
Gd in the diquark channel. The diquark mass Md(0, 0) in the vacuum is determined 

by 


1-8G 


d^p 


1 


1 


° J {27T)3\Ep + MDiO,0)/2^ Ep-Md{0,0)/2) 

In general, temperature effect disorders a system, and any condensate will be 
suppressed in hot medium. Therefore, the critical chemical potential of color super¬ 
conductivity goes up monotonously with increasing temperature, and the minimum 
one /iA at zero temperature satisfies the gap equation 


1-8Gd 


d^p / 1 

(27r)3 VEp -I- 



( 10 ) 


From the comparison of the above two equations, we obtain the relation between 
the diquark mass in the vacuum and the critical chemical potential for diquark 
condensate at zero temperature, /ta = M£)(0,0)/2 . 

When the diquark is massless in the vacuum, the color superconductivity hap¬ 
pens even at ^a = 0, and the vacuum becomes instable. On the other hand, when 
the diquark mass in the vacuum is larger than two times the quark mass, the color 
superconductivity is far from the vacuum, and a diquark can not be in the bound 
state, but decay into two quarks. Therefore, the diquark mass Md{0, 0) should sat¬ 
isfy the constraint 0 < Md{0, 0) < 2Mq( 0, 0) . From the comparison of the diquark 
mass equation m and the quark mass equation in the vacuum. 


1 - 24Gs 


d^p I 

(2^)" ^M2(0,0) + p2 


( 11 ) 


the constraint 0 < M£)(0,0) < 2Mq{0, 0) on the diquark mass resulus in the relation 
between the two coupling constants Go and Gs, 


G min ^ ^ r^max 

D ^ 5 


4 (^A^M^(0,0) + + M^(0,0) In ) ’ 

^3 Mq(0,0) 

2 ^MQ(0,0)-m’ 


( 12 ) 


With the known values of the parameters A,to and Mq( 0, 0) in real world, the 
minimum and maximum coupling constant take the values GJ)*" = 0.82Gs and 

G-rnax ^ 1 55 ^^ 


4. Phase Diagram 

Before making numerical calculations, we first discuss the parameters in the model. 
There are four independent parameters in the current NJL model with both diquark 
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and meson channels, the three-momentum cutoff A, the current quark mass to, and 
the coupling constants Gs and Gd- The first three can be fixed by fitting the pion 
mass TO,r = 0.134 GeV, the pion decay constant /,r = 0.093 GeV, and the quark 
condensate density a = (—0.25 GeV)^ in the vacuum. It leads to A = 0.65 GeV 
and Gs = 5.01 GeV“^ in chiral limit with zero current quark mass to = 0, and 
A = 0.653 GeV, Gs = 4.93 GeV“^ and to = 0.005 GeV in real world. 

In chiral limit the chiral phase transition is well dehned, and the phase transition 
line determined by cr = 0 is a solution of the gap equations The phase diagram 
in the T—ft plane is shown in Fig^for different coupling constant in the diquark 
channel. 

At the minimum coupling constant G^*”, the whole T—ft plane is separated into 
three regions, the region with chiral symmetry breaking but color symmetry, cr 7 ^ 0 
and A = 0 at low density, the region with chiral symmetry but color symmetry 
breaking, cr = 0 and A 7 ^ 0 at low temperature, and the region with both chiral and 
color symmetries, ct = A = 0 at high temperature and/or high density. The critical 
temperature = 185 Me V for ch iral phase transition at /i = 0 is in good agreement 
with the lattice calculatior^El and the chiral symmetry breaking phase ends at 
fi^ determined by 

l-24Gs/i0rh(t-^,) = O. (13) 

From the comparison with the quark mass equation lllll in the vacuum, one has 
fty- = Mq( 0, 0). Therefore, from the dehnition of the minimum coupling con¬ 
stant G^®” at which the minimum chemical potential ft^ of color superconduc¬ 
tivity is equal to the constituent quark mass in the vacuum, the end point of chi¬ 
ral phase transition and the starting point of color superconductivity coincide at 
p = Mq( 0 , 0 ). 

The order of a phase transition is determined by th e b ehavior of the order 
parameter across the phase transition line. It is well knowrP^that the chiral phase 
transition is of second order at high temperature and of first order at high density. 
The critical chemical potential at which the chiral transition goes from first to 
second order is 285 MeV in the NJL mode^^. One can also determine the critical 
exponent for the color superconductivity. Since the square bracket of the second 
gap equation of 0 depends quadratically on the order parameter A, one has a 
simple zero for A^ as a function of T and therefore in the neighborhood of the 
critical temperature Ta of color superconductivity the order parameter behaviors 
as A(r) ~ |T — TaI^^^ . This means that the transition between the phase with 
(T = 0, A 7^ 0 and the phase with cr = A = 0 is of second order. 

With increasing coupling constant Gd, the critical temperature Ta for melting 
the diquark condensate increases and the critical chemical potential /iA for forming 
the color superconductivity decreases. Especially, a new phase with both chiral and 
color symmetry breaking, cr 7 ^ 0 and A 7 ^ 0, is produced in between the phase with 
only chiral symmetry breaking and the phase with only color symmetry breaking. 
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as shown in Fig.?? for Gd = 1.4Gs. The left and right borders of this mixed phase 
are determined by ct 0, A = 0 and ct = 0, A ^ 0, respectively. With increasing 
Gd the chiral breaking phase is swallowed up gradually by the mixed phase. At the 
maximum coupling constant the chiral breaking phase disappears and the 

mixed phase reaches the maximum. 





Fig. 1. The Phase diagram in the T — (i plane for different coupling constant Gd ihe diquark 
channel. 


5. Summary 

We have analytically derived the simple relation between the diquark mass in 
the vacuum and the minimum quark chemical potential for diquark condensate, 
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fiA = M£){0,Q)/2. Under the conditions to keep the diquark bound state in the 
vacuum and to have a stable physical vacuum without diquark condensate in 
it, we obtained the constraint on the coupling constant in the diquark channel, 
Qmm Qjj G^““, with the relations to the coupling constant in the meson chan¬ 
nel — 0.82Gs and G2)““ — 3/2Gs. There is a strong competition between the 

color superconductivity and chiral symmetry restoration: With increasing Gd the 
phase with chiral symmetry breaking but color symmetry is gradually taken over 
by the phase with both chiral and color symmetry breaking. 
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